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We investigate effects of short range interactions on the addition spectra of quantum dots using 
a disordered Hubbard model. A correlation function S{q) is defined on the inverse compressibility 
versus filling data, and computed numerically for small lattices. Two regimes of interaction strength 
are identified: the even/odd fiuctuations regime typical of Fermi liquid ground states, and a regime 
of structureless S[q) at strong interactions. We propose to understand the latter regime in terms of 
magnetically correlated localized spins. 
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Coulomb interactions and disorder in electronic sys- 
tems have posed a major challenge to condensed matter 
physics for quite some time. 

Quantum dots with discrete electronic spectra offer a new 
avenue to this problem. A direct probe to the ground 
state energy is given by Coulomb blockade peaks in the 
conductance as the gate voltage is varied Theory 
of spectral fluctuations of non interacting electrons has 
made much progress during the last decade due to the 
advent of semiclassical approximations, random matrix 
theory and the non linear sigma model approach 
However since Coulomb interactions are essential for the 
"Coulomb blockade" effect, one may wonder as to the 
validity of non interacting approximations to quantum 
dots in general. In particular: Is the ground state quali- 
tatively similar or different than a Fock state of the lowest 
single electron orbitals? 

To gain insight into this question, we consider a system 
of interacting electrons on a finite tight binding lattice 
with onsite disorder. The inverse compressibility at con- 
secutive fillings is 



A(A^) = E{N + 1) - 2E{N) + E{N - 1) 



(1) 



where E{N) is the ground state energy of a dot with N 
electrons. (We assume weakly coupled leads such that 
is well defined within the area of the dot.) By vary- 
ing a gate potential t^, the dot's energy is modified to 
E"^ = E{N) — eifN. Conductance peaks through the 
leads are observed at E'^iN) = E'^iN + 1), i.e. at po- 
tentials etpN = E{N + 1) - E{N). Thus differences be- 
tween the peak potentials (pN yield direct measurements 
of A(A^) which can be defined as times the discrete 
inverse capacitance of the dot. 

We shall model the single electron part of the dot's 
Hamiltonian by a site-disordered tight binding model 



(ij) 



E 



(2) 
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where c^^ creates an electron at site i with spin s, (ij) 
denote nearest neighbors on the lattice, and Wi are ran- 
dom site energies taken from a uniform distribution in 
the domain [— M^/\/3, Vl^/\/3]. aJis creates an electron 
in eigenstate 0„ and spin s. An orbital coupling to a 
magnetic field is included by defining tij {B) = te^'^^-^ii ^ 
where V x A = i?. 

In the absence of electron interactions, the inverse com- 
pressibility is given by 



Ao(Ar) = 







A^ = 2n -I- 1 



e„+i - e„ N = 2n 



(4) 



where e„ are defined in (^). We find it useful to define 
the "Ising" variables 



S{N) 



A{N) - A{N - 1) 
|A(Ar)-A(A^-l)| 



(5) 



and a corresponding correlation function on a series of L 
consecutive data points, 

1 ^ 

'5(9) ^ 72 E S{N..)S{N,) exp (-*(A^, - N,)q) (6) 



i,3 = l 



Obviously, the non interacting spectrum has perfect 
"long range antiferromagnetic correlations" i.e. iSo(7r) = 
1. 

Coulomb interactions are treated by separating the in- 
teractions into the long and short range parts. A crude 
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approximation to these two terms is given by an onsite 
and an infinite range parts 



Hi, 



2C 



(7) 



where Uis = c^gCi^, s =Ti i- It is clear that the first term 
simply adds a constant ^ to A(A^), and therefore does 
not alter S{q). Thus in our model, deviations of S{q) 
from Sq {q) must therefore be a consequence of the Hub- 
bard interactions described by the second term in (^. 

We restrict ourselves to a square lattice of Af sites with 
periodic boundary conditions, and to disorder strength 
W appropriate for the "diffusive" regime, i.e. the mean 
free path I for the non interacting electrons is of the order 
of, or smaller than, the system's linear length L. Using 
I = vpT, where the inverse lifetime is calculated in the 
Born approximation to be t = r>^w'^ ' ^® ^'^"^ L/l ^ 

Perturbation theory: We diagonalize Ti.Q on a square lat- 
tice of size Af, with periodic boundary conditions for a 
given realization {wi}. The single electron eigenenergies 
{e„} and wavefunctions ipn{i) are assumed to be known. 
The first order correction to Eq. (^) are given by second 
differences of the first order energies 



(8) 



Here we appeal to the random matrix properties of Ti.Q 
in the diffusive regime, in order to estimate the magni- 
tude and fluctuations of Ei analytically. We assume a 
Random Vector Model (RVM) where all eigenvectors ipn 
are random complex unit vectors of dimension TV whose 
ensemble averaged correlations are [Q 



r 



(9) 



Using the orthonormalization constraints for ipn we ob- 
tain after some algebra |^ that the RVM estimate for the 
average first order correction to A is 



A™™ (AT) 



N odd 



3U 
71+2 

-M2 ^even 



(10) 



We have compared the RVM estimates to numerical re- 
sults for disorder averaged A""™(iV) for odd and even 
N respectively. Calculations have been done for lattice 
sizes Af ~ 56,110,210,420, with the disorder varied in 
the range l/L e [0.1,2.5]. We find that 



A7™(7V) e 



[2.3,3.5] X ^ 
[-1.3,-2.71 X 



N odd 
N even 



(11) 



Comparison of (10) to ( [ll[ ) shows the RVM estimates to 
be in the right ball park as those gotten by numerically 



determining ipn of the disordered tight binding model. 
The main lesson learned by this calculation is that first 
order Hubbard corrections reduce on average the fluctu- 
ations in A{N) since they are positive for odd N and 
negative (on average) for even N. 

Variational Theory for Sign Flips. The weak coupling 
regime is defined where the non interacting ground state 
1 4*0) is variationally stable against particle hole exci- 
tations. When interaction strength exceeds a certain 
threshold, it is variationally advantageous to create spin 
polarized electron-hole pairs e.g. c]^p^_|_;^c„^ J^o) which 
reduce the Hubbard interaction energy at the expense 
of enhanced single particle (kinetic) energy e„p+i — Cnp- 
In this variational theory, for an even number of elec- 
trons the threshold for forming a triplet is given by the 
inequality 

> Ae - {El{2n) - El{2n)) ^ Ae + UT 

:F^J2\ i\^n+l{^^)\' - |^„(X.)|') (j2 IV-n'l^On 



(12) 



where Ae = e„+i — e„ and E\,El are the interaction 
corrections to the singlet and triplet energies respectively. 
Eq. ( p^ ) can alternatively be written as an inequality for 
A(A^) which includes up to first order corrections in U : 

0> A(2n) - A(27i - 1) 

- (|^„+i(x,)|2 - |V'„(x,)|2) |7A„(x,)|' (13) 



For extended random wave functions, the last term in 
( p^ is readily seen to be of order Af~^^'^ and thus negli- 
gible in comparison to A{2n) — A(2n — 1). Eqs. (12 1^) 
establish the connection between formation of triplets 
and sign flips of SiN) (defined in (|)). The sign flips 
degrade the "antiferromagnetic" correlations of S{q) as 
U is increased. We shall proceed to estimate the leading 
dependence of S{tt; U) by a statistical calculation based 
on the properties of the non interacting spectrum of Hq 
with and without an orbital magnetic field. 

We assume that the level spacing statistics of Hq in the 
diffusive regime of L/Z ~ 1 is that of a random matrix in 
the Gaussian Orthogonal Ensemble (GOE). In the pres- 
ence of a magnetic field B whose flux through the dot is of 
order of one flux quantum divided by Af^^'^ ^ , the level 
spacing statistics of Tio turns into the Gaussian Unitary 
Ensemble (CUE). These assumptions were checked nu- 
merically and verified quite well for the Hamiltonian (|^) 
with and without external flux, on lattices up to A/" = 420 
sites. 



The Wigner distributions of Ae are 0] 
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(Ae) = 
(Ae) = 



2A2 



32(Ae)2 4(A.)^ 

= — e tta^ 

7r2A3 



(14) 



where A is the mean level spacing. 

The probablity distribution of the interaction term in 
( p!^ is denoted by 



(15) 



where P{ is a dimensionless function of its dimensionless 
argument. The probabihty of satisfying Eq. (|l^) is given 
by the double integral 



pfiw^jj) ^ / dTP[{T) / d(Ae)P2(Ae) (16) 
Jo Jo 



which at weak coupling U << A we can evaulate using 
the low energy expansion of P2 and obtain 



pL% - J- d.P[{.) /; dee = U^~{^') (17) 



4A2 
32 



for = and 5(7r) = 1 — {).2U^. This supports the va- 
lidity of the GOE (GUE) statistics for the case (f) — Q 
(<^ = V2). 

Experimental S{q). In Fig. (0) we have plotted three 
sets of data for A(iV) and S{q). The top two sets were 
measured on two GaAs samples at different fillings. The 
bottom data, measured on InO, was taken from Ref. pO[ . 
It is interesting to note that S{q) shows strong even- 
odd correlations only for the middle data set, while it 
has no remnance of these correlations for the top and 
bottom sets. Although we do not attempt to provide a 
quantitative understanding of the data in this paper, we 
dare to speculate that the flatness of S{q) correlations is 
definitevely associated with strongly correlated samples 
where ground states significantly differ from Fermi-liquid 
Fock states. The fluctuations in A(A^) due to the direct 
Coulomb interaction and its interplay with disorder are 
discussed in Ref. [ pT] . 
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At weak coupling the reduction in the perfect even-odd 
correlations are proportional to P^'^p, and therefore we 
find that 
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5(7r) cx 
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GOE 
GUE 



(18) 



Numerical Diagonalizations. The variational theory and 
the random matrix estimates were checked by numeri- 
cally diagonalizing the full interacting Hamiltonian on a 
lattice size 3x3 with periodic boundary conditions with 
a magnetic flux (p threading the lattice at its center. We 
varied the number of electrons for each specific realiza- 
tion from = 5 to TV = 13, each time diagonalizing the 
corresponding (^) matrix (i.e., maximum matrices of 
size 48620 x 48620 for TV = 9) for various values of U. 
Then ^(Tr) was calculated and averaged over 500 different 
realizations. 

Two different values of (p were considered: ((> — and 
(j) = tt/2 which should correspond to the GOE and GUE 
statistics. Because of the small size of the samples the 
maximum flux oi (j) = 7r/2 is not strong enough to com- 
pletely remove time reversl symetry, and the level spacing 
is a combination of the GOE and GUE level spacing. The 
results of S{Tr) for both cases are presented in Fig. ||. It 
can be seen that the general behavior predicted in Eqs. 
( p^ ) is observed. A reasonable fit for the numerical re- 
sults in the regime C/ < 1 is obtained by S{tt) — 1 — 0.3U^ 
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FIG. 1. Averaged S{q), Eq. for the disordered Hub- 
bard model on a 3 X 3 torus. Inset: S{tt) for the (p = and 
= 7r/2 cases, fitted to estimates of RVM variational theory. 



FIG. 2. Experimental second energy differences from quan- 
tum dots of GaAs, Ref. Q (top and middle), and InO, Ref. 
[0, (bottom), with corresponding correlations S{q). Corre- 
lations indicates that the top and bottom datae are in the 
strongly interacting regime. 
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